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In the presence of external, linear / nonlinear electromagnetic fields we integrate f{R) ~ -R + 
2a%/ R + const, gravity equations. In contrast to their Einsteinian cousins the obtained black holes 
are non-asymptotically flat with a deficit angle. In proper limits we obtain from our general solution 
the global monopole solution in f{R) gravity. The scale symmetry breaking term adopted as the 
nonlinear electromagnetic source adjusts the sign of the mass of the resulting black hole to be 
physical. 
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' / {R) gravity is a modified version of standard Einstein's gravity which incorporates an arbitrary function of the 
h-H ' Ricci scalar (i?) instead of the linear one [l| (for a recent review). Depending only on the Ricci scalar may sound 
\Q I simpler at the initial but the pertinent nonlinearity makes nothing simpler than the Einstein's gravity with sources. 
There are both advantages and disadvantages in adopting such a model. It contains for instance, its own source 
known as the curvature source in the absence of an external matter source. The identification of physical sources, 
^ ' however, within the nonlinear structure through its equations is not an easy task at all. For the same reason almost 
all known solutions, except very few, result in nonanalytical (i.e. numerical) expressions for the function / (i?) . 
^ Starting from a known function of / (i?) a priori is an alternative approach which hosts its own shortcoming from 
j^JT), the outset. Keeping a set of free parameters to be fixed by observational data can be employed in favour of / [R) 
'— gravity to explain a number of cosmological phenomena. First of all, to be on the safe side along with the successes of 
' general relativity most researchers prefer an ansatz of the form / (i?) = R + ag (R) , so that with a — > one recovers 
^ the Einstein limit. The struggle now is for the new function g (i?) whose equations are not easier than those satisfied 
by / (i?) itself. Without seeking resort to this latter (and easier) route we have shown recently that / (R) = ^/R 
OO I gravity admits exact solution in 6— dimensional spacetime with the external Yang-Mills field 0. Without demanding 
' an analytical representation for / (i?) , as a matter of fact, exact solutions are available in all dimensions with the 
Yang-Mills source. Similar results may be investigated with other sources such as the Maxwell fields. This will be 
our strategy in the present Letter. 
_ We assume / (i?) = {R + Ri ) + 2a^/R + Rq , in which ^, a, Rq and Ri are constants, a priori to secure the_Einstein 
• limit by setting the constants Rq = Ri = a = and ^ = 1. This extends a previous study without sources [sj to the 
I ' case with sources. Why the square root term in the Lagrangian?. It will be shown that for i?o = = and without 
^ . external sources such a choice of square root Lagrangian gives the curvature energy-momentum tensor components 
as Tj* = , Tg = — 0, which signify a global monopole '4*|. A global monopole which arises from spontaneous 
^ . , breaking of gauge symmetry is the minimal structure that yields non-zero curvature even with zero mass. We test 
' the analogous concept in f{R) gravity to obtain similar structures. Unlike the case of our concern here will be 
restricted to the 4— dimensional spacetime. As source, we take electromagnetic fields, both from the linear (Maxwell) 
and the nonlinear theories. For the linear Maxwell source we obtain a black hole solution with electric charge (Q) 
and magnetic charge (P) reminiscent of the Reissner-Nordstrom (RN) solution with different asymptotic behaviors. 
That is, our spacetime is non-asymptotically flat with a deficit angle. For the nonlinear, pure electric source we 
choose the standard Maxwell invariant superposed with the square root invariant, i.e. the Lagrangian is given by 
£ {F) ^ F + 2/?V— P, where F = jF^^^F^'^ is the Maxwell invariant and /3 is a coupling constant. This particular 
choice has the feature that it breaks the scale invariance [5| , gives a linear electric potential which plays role in quark 
confinement [1, We find out that the scale breaking parameter (3 modifies the mass of the black hole. For this 
reason Lagrangians supplemented by a square-root Maxwell Lagrangian may find rooms of applications in black hole 
physics. 
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II. / {R) GRAVITY COUPLED WITH MAXWELL FIELD 



The action for / (i?) gravity coupled with Maxwell field in 4-dimensions is given by 

7(^) 1 
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s = / <rx4^ 



2k 47r 
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in which / (i?) is a real function of Ricci scalar R and F = jF^j^F'^'^ is the Maxwell invariant. (We choose k = Sw 
and G = 1). The Maxwell two- form is chosen to be 



Q 



dt A dr + P sin ede A 



(2) 



in which Q and P are the electric and magnetic charges, respectively. Our static spherically symmetric metric ansatz 
is 



dr 



ds^ = -A (r) dt^ + + {dO^ + sin^ ed(tf) 



(3) 



where A (r) stands for the only metric function to be found. The Maxwell equations (i.e. dF = = d*F) are satisfied 
and the field equations are given by 



in which 



dR ' 
1 



while the energy momentum tensor is 



47rT; = -FS'; + F^xF"^. 



Furthermore, the trace of the field equation (4) reads 
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fRR+{d-l)nfR--f = kT 
with T = Tj^. The non-zero energy momentum tensor components are 

p2 + g2 . 
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with zero trace and consequently 



One finds 



f = ^fRR + SDfR. 
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in which a prime denotes derivative with respect to r. Overall, the field equations read now 
Herein 

and for the details we refer to Q. The tt and rr components of the field equations imply 
or equivalently 

/^ = 0. 

This leads to the solution 

where ^ and 77 are two positive constants Q . The other field equations become 



Now, we make the choice 



which leads to 



1 



R — ^2^2 ^ 



where a, Rq and ^ (from (21)) are constants. As a result one obtains for / (r) 



/ - H TTt-Ko 



and from (12) we have 

r'^A" + 4rA' + 2 (yl - 1) 



Rq. 



This equation admits a solution for the metric function given by 

2?7^ r r-^ 12 

Herein the two integration constants Ci and C2 are identified through the other field equations (22) and (23) 
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while for the free parameters we have a = t] > 0. Finally the metric function becomes 

A{r)^--- + \-^r^ 30 
2 r 3 

where m = — < 0, Ae// — —j^ and = - — | — -. The choice of the free parameters in terms of each other 
prevents us from obtaining the general relativity limit, namely the Reissner-Nordstrom (RN)-de Sitter (dS) solution. 
It is observed that the parameter ^ acts as a scale factor for mass and charge and for the case £, = 1 and Q = P = the 
solution reduces to the known solution given by 0, • The properties of this solution can be summarized as follow: 
The mass term has the opposite sign to that of Schwarzschild and the solution is not asymptotically flat, giving rise 
to a deficit angle. The latter property is reminiscent of a global monopole term with a fixed charge. To see the case 
of a global monopole we set Ra — Q = (i.e. zero external charges and zero cosmological constant) and find the 
energy-momentum components. This reveals that the non-zero components are Tf = = which identifies a 

global monopole [1]. The solution (30) can therefore be interpreted as an Einstein-Maxwell plus a global monopole 
solution in f{R) gravity. The area of a sphere of radius r (for q^ = Rq = 0) is not Airr^ but 27rr^. Further, it can 
be shown easily that the surface ^ = f has the geometry of a cone with a deficit angle A = -I [3|a. It can also be 
anticipated that a global monopole modifies perihelion of circular orbits, light bending and other physical properties. 
Although in the linear Maxwell theory the sign of mass is opposite, in the next section we shall show that this can be 
overcome by going to the nonlinear electrodynamics with a square root Lagrangian. Another aspect of the solution is 
that since f^. > we have no ghost states. 



III. f{R) GRAVITY COUPLED WITH NONLINEAR ELECTROMAGNETISM 
A. Solution within nonlinear electrodynamics 

In this section we use an extended model for the Maxwell Lagrangian given in the action 

7(i?) 
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where / (i?) = ^ (i? + i?i) + 2ay/R + Rq, in which i?i and Rq are constants to be found while 

CiF)^-l-{F + 2pV^). 



(31) 



(32) 



Here /3 is a free parameter such that lim^_j.o C (F) — ~j^F, which is the linear Maxwell Lagrangian. The main reason 
for adding this term is to break the scale invariance and create a mass term Q . The normal Maxwell action is known 
to be invariant under the scale transformation, x — >■ Xx, A^^ — ^ x^m' (-^ =const.), while ^/—F violates this rule. We 
shall show how a similar term modifies the mass term in f{R) gravity. Our choice of the Maxwell 2-form is written as 



F ^ E{r)dtA dr 

and the spherical line element as (3). The nonlinear Maxwell equation reads 



d *Y- 



dF 







which yields the solution 



E{r)^^/2p + % 



(33) 
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with a confining electric potential as V {r) — — \/2Br + ^ . This is known as the " Cornell potential" for quark 
confinement in quantum chromodynamics (QCD) [^Q- The Einstein equations implies the same equations as (4-7) 
and the energy momentum tensor 



t; = c{F)s;-F^xF 



dC 

'of 



(36) 
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with the additional condition that the trace = T 7^ 0, here. Upon substitution into the field equations one gets 

Ri = ^+2^0. (37) 
a = r] (38) 
and a black hole solution results with the metric function 



"^^'^=2- Srjr +e^ + T2"- ^^^^ 



This is equivalent to the solution given in (30) with the same ^eff but with the new m = '^^^^ ^ and q — This 
is how the scale breaking term in the Lagrangian modifies the mass. 

For the sake of completeness we comment here that, choosing a magnetic ansatz for the field two- form as 

F = PsinOde Adip (40) 

together with a nonlinear Maxwell Lagrangian 

^(^)--^(^ + 2/3Vf) (41) 



and 



Ri = ^Ro (42) 



admits the magnetic version of the solution as 

1 472/3P-e Ro^2 

2 3rjr ^r'^ 12 

The magnetic solution, however, is not as interesting as the electric one 



^ (^) - o - o^. + 712 + TT^'- (43) 



B. Thermodynamical aspects 

The solution we found in the previous section is feasible as far as a physical solution is concerned. Here we set our 
parameters, including the condition ^ and r/ positive, to get 4-\/2/3(5 — £, > such that the solution admits a black 
hole solution with positive mass as 

^i^)= + ^ + 44 

2 r 12 

Now we wish to discuss some of the thermodynamical properties by using the Misner-Sharp % 1^ energy to show 
that the first law of thermodynamics is satisfied. To do so first we set Rg — and introduce the possible event horizon 
as r = r/i such that A{rh) = 0. This yields 



in which 



r± = m± \J m2 - 2q^ (45) 
^Th = f+) 

A(r)^(^-^:4$-^ (46) 
2r^ 

and the constraint m > m^ri is imposed with merit — V^q ■ If one sets Q > 0, this condition is satisfied if 
Q > ——> — ^ — (providing 4/3 + — |- ^ 0). The choice m = merit leads to the extremal black hole. The Hawking 

temperature is defined as 



A' (r+) _ rl -^2^ 

.3 
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and the entropy 



^ fnl 



(48) 



with Aj, 

by 



47rr^, the surface area of the black hole at the horizon. The heat capacity of the black hole also is given 
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(49) 



which takes both (+) and (— ) values. Both the vanishing / diverging Cq values indicate special points at which the 
system attains thermodynamical phrase changes. 
The first law of thermodynamics can be written as 



TdS -dE = PdV 



in which 



dE = — 
2k 
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if - Rfu) 
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— /fl; stands for the black hole 



with E the Misner-Sharp energy and T = ^ the Hawking temperature. Further, S — ^ 
entropy, p = = Tq is the radial pressure of matter fields at the horizon and finally the change of volume of the 
black hole at the horizon is given by dV — A+dr-^-. One can easily show that the first law in thermodynamics in the 
form introduced above is satisfied. 



IV. CONCLUSION 



Exact solutions for nowadays popular, modified gravity model known as f{R) gravity with external sources (i.e. 
Tpf-**^^ 7^ 0) are rare in the literature. We attempt to fill this vacuum partially by considering external electromagnetic 
fields (both linear and nonlinear) in /(i?) gravity with the ansatz / (i?) = ^ (i? + i?i) + 2a\/R + Rq. In this choice Rq 
is a constant related to the cosmological constant, the constant Ri is related to i?o while a is the coupling constant 
for the correction term. This covers both the cases of linear Maxwell and a special case of power-law nonlinear 
electromagnetism. The non-asymptotically fiat black hole solution obtained for the Maxwell source is naturally 
different and has no limit of the RN black hole solution. In the limit of Q = P = Ag/ / = we obtain the metric 
for a global monopole in f{R) gravity. Our solution can appropriately be interpreted as a global monopole solution 
in the presence of the electromagnetic fields. The thermodynamical properties of our black hole solution is analyzed 
by making use of the Misner-Sharp formalism and shown to obey the first law. As the nonlinear electromagnetic 
Lagrangian we choose the normal Maxwell, supplemented with the square root Maxwell invariant which amounts to a 
linear electric field. This latter form is known to break the scale invariance yielding a linear potential which is believed 
to play role in quark confinement problem. Within f{R) gravity the presence of scale breaking term modifies the mass 
of the resulting black hole. The advantage of employing square-root Maxwell Lagrangian as a nonlinear correction 
can be stated as follows: Beside confinement in the linear Maxwell case we have in f{R) gravity an opposite mass 
term while with the coupling of the aquare-root Maxwell Lagrangian we can rectify the sign of this term. 
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